Magnetization switching by ultrashort acoustic pulses 
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It is shown theoretically that a single a lew picoseconds long acoustic pulse can reverse mag- 
netization in a magneto-strictive material Terlenol-D. Following giant magneto-elastic changes ol 
free energy density the magnetization vector is kicked out ol a local in-plane energy mininum and 
decays into another minimum. For acoustic pulse duration significantly shorter than magnetization 
precession period r ac <C T prsc , the switching threshold is determined by the acoustic pulse area, 
i.e. pulse integral in time domain, similar to coherent phenomena in optics. Simulation results are 
summarized in a magneto-acoustic switching diagram and discussed in the context ol all-optical 
magnetization switching by circularly polarized light pulses. 
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Searching for new possibilities of ultrafast magnetiza- 
tion switching is motivated by continuously growing de- 
mand in faster data recording technologies, which are 
based on reversal of individual magnetic bits at the nano- 
scale. Among different mechanisms of magnetization 
switching the acoustically driven switching at ultrafast 
time scales remains largely unexplored. 

The first pioneering time-resolved observation of mag- 
netization dynamics in ferromagnetic nickel induced by 
femtosecond laser pulses revealed the phenomenon of 
ultrafast demagnetization in nickel The ultrafast 

drop of magnetization on a sub-picosecond time scale 
was caused by transient elevation of electron tempera- 
ture close to the Curie point. The subsequent dynam- 
ics of magnetization recovery on a time scale exceed- 
ing tens of picoseconds could be adequately reproduced 
by Landau-Lifshitz-Gilbert (LLG) equations describ- 
ing damped precession of magnetization vector in the 
presence of temperature-dependent magneto-crystalline 
anisotropy Q. 

The most recent experiments combining picosecond 
acoustics [H with ultrafast magneto-optics showed that 
the magneto-crystalline anisotropy can be also changed 
by picosecond acoustic strain pulses to a ferromagnetic 
semiconductor GaMnAs [3, |B| or nickel 0, thus trig- 
gering the magnetization precession without heating the 
sample. However, due to the relatively small magneto- 
strictive coefficient in GaMnAs (A ~ 10~ 5 at cryogenic 
temperature [7]) and nickel (A - 3 x 1(T 5 at room tem- 
perature [![), the resulting transient change in magneti- 
zation direction appeared to be too small and the magne- 
tization vector returned to its initial direction, i.e. mag- 
netization switching (reversal) did not occur. 

The first experimental demonstration of ultrafast non- 
thermal magnetization switching was reported in 1998 
by Back and co-workers Q, who used ultrashort pulses 
of magnetic field induced by relativistic electron bunches 
to switch magnetization in Co/Pt film. Most recently 
a spectacular observation of all-optical magnetization 



switching in GdFeCo using circularly polarized light 
pulses O, El, raised a lot of questions about physically 
unclear switching mechanisms suggesting that not only 
transient overheating of electrons but also the dynamics 
of lattice temperature may be necessary to explain the 
underlying physics 1 1 la . Not only the minimum amount 
of deposited heat [12[ but also a minimum amount of 'cir- 
cularity' was necessary to obtain switching 11311 . at least 
within a certain range of optical excitation [14j . 

In this letter we theoretically investigate the in- 
teraction of ultrashort acoustic pulses with Terfcnol- 
D (Tb x Dyi_ 2; )Fc2, the rare-earth compound famous 
through its giant magneto-strictive coefficient A ~ 10~ 3 
8] , and demonstrate the possibility of ultrafast magneto- 
acoustic switching. Moreover, the results are discussed 
within the framework of recently observed all-optical 
magnetization switching by single circularly polarized 
femtosecond laser pulses, which is undoubtedly accompa- 
nied by the generation of picosecond pulses of coherent 
acoustic phonons Q. 

The phenomenological expression for free energy den- 
sity F$) = F k + F me + F d + F z for (110) thin films of 
Tbo.27Dyo.73Fe2 epitaxially grown on sapphire reads fljj ]: 



F k = 

TP — 
1 m.p. — 



F d = 



K\{p? x a 
bi{a 2 x e xi 
b2(a x a y e xy 

^(M s cos0 



2 2 

OZ.O& 



a y e 



yy 
a x a : 



2 2\ 

a z a y ) 
- die. 



K 2 {ala 2 y al) 



(1) 

(2) 
(3) 



where a x , a y , a z are direction cosines of the magnetiza- 
tion vector M = M s (a x ) in the crystallographic 
coordinate system (x, y, z) and saturation magnetiza- 
tion HoM s — 0.945 T. In the rotated frame (x',y',z') 
the direction of magnetization is determined by two 
angles 9 (out-of-plane angle) and (f> (in-plane angle), 
see Fig. 1(a). In Eq. (fl][3j) F k , F me and F d denote 
the magneto-crystalline anisotropy, magneto-elastic and 
magneto-static terms, respectively, and Zeeman contri- 
bution F z = — /io^ea;*^ is disregarded through out this 
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paper as we consider the case of zero external magnetic 
field, ~fi ext = 0. 

Epitaxial growth of a thin Terfenol-D film in (110) 
direction on a lattice-mismatched substrate induces the 
built-in static strain described by the following tensor: 
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which is determined by a single strain component e xy —- 
0.55% in a crystallographic coordinate frame {x, y, z) 
(Fig. 1(a)). The competition of different contributions 
in the total free energy density results into four local in- 
plane energy minima corresponding to four different mag- 
netization directions 1,2,3,4 in Fig. 1(b). The explicit 
dependence of the magneto-elastic term F me both on 
the strain components and magnetization direction and 
large values of magneto-elastic coupling coefficients b\ = 
— |A 100 (cii_— ci 2 ) = —80 MPa and b 2 = — 3A ni c 4 4 = 



—85 MPa 17] suggest that application of external strain 
will shift the minima of free energy minima and there- 
fore, change the magnetization direction. Indeed, the 
application of time-independent uniaxial strain 77 in the 
direction normal to the surface of a thin Terfenol-D film 
results into the in-plane shift of all four energy minima 
by the angle A$, as illustrated in Fig. 1(c) for rj — —0.3% 
(film compression) and r\ = 0.9% (film tension). 

The action of a time-dependent uniaxial strain (acous- 
tic pulse) r] (t) can be described by adding the following 
dynamic tensor in a crystallographic coordinate frame 
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where the rotation of coordinates system by 45° from 
(x',y',z') into (x,y,z) leads to factor 1/2 and generates 
the non-diagonal terms. 

Inserting the total strain e(t) = e s tat + ed yn (t) in 
Eq. ((2]) generates explicit time-dependence of free energy 
F(t), which drives the magneto-acoustic dynamics de- 
scribed by LLG equation [6| 

- = - 7 ^x^ //) + -^x-j (6) 

where the first term describes the torque driving the pre- 
cession of magnetization vector around the effective time- 
dependent magnetic field H e ff(t) 
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and the second term describes precession damping ac- 
cording to the phenomenological Gilbert damping param- 
eter a = 0.1 [18|; 7 is the gyromagnetic ratio. 
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FIG. 1: Magneto-elastic free energy density in a thin (110) 
film of Terfenol-D grown on sapphire (a) possesses four in- 
plane energy minima (b). Application of a static uniaxial 
compressive (77 = —0.3%) or tensile (77 = 0.9%) strain in the 
direction perpendicular to the film leads to the in-plane shift 
A<E> of all four energy minima (c). 



If the system is initially prepared in minimum 2, the 
application of instantaneous step-like strain [l9[ results 
into damped precession of magnetization vector into the 
new minimum 2' with a precession period T prec ~ 25 ps, 
see Fig. 2. However, if the strain is turned off after some 
time, the situation equivalent to the application of a rect- 
angular acoustic pulse of finite duration r ac , the preces- 
sion trajectory will decay back into minimum 2. Such 
magneto-acoustic precession trajectory induced by the 
action of a picosecond acoustic pulse with r ac = 3 ps and 
strain amplitude r/ ac — 3 x 10~ 3 is shown in Fig. 2 and 
can be explained analytically. 

LLG equations (J5]) for the initial condition in one of 
the four energy minima read 
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FIG. 2: Starting from energy minimum 2, a dynamic step- 
like strain induces a damped magnetization precession around 
the shifted minimum 2' (dashed line). An ultrashort acous- 
tic pulse results into the out-of-plain kick A0 ac followed by 
damped precession into the initial minimum 2 (continuous 
line). 



and their integration for an ultrashort acoustic pulse rj(t) 
obeying the condition r ac <C T prec approximate well the 
acoustic out-of-plane kick of the magnetization vector by 
the angle 
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The product of directional cosines ct' y a' z in Eq. (JSJ equals 
0.48 for energy minima 1 and 3 and -0.48 for two other 
minima, respectively. Therefore, depending on the initial 
condition, the same acoustic pulse will kick the magneti- 
zation vector out of sample plane in opposite directions. 
Equation Q clarifies the microscopic physical model be- 
yond the time-dependent magnetic torque 1 5? x H e ff | in- 
troduced by Kim and co-workers @ and shows that the 
prefactor in Eq. ([9]) is dominated by the ratio of magneto- 
elastic coupling coefficients &i and 62 (which are both 
proportional to the respective coefficients of magneto- 
strictive tensor A) to saturation magnetization M s . 

Similar to polarization dynamics in coherent optics, 
the acoustic rotation angle A8 ac of magnetization vec- 
tor appears to be proportional to the acoustic pulse 
area J rj(t)dt for arbitrary acoustic pulses obeying r ac <C 
T prec ~ 25 ps. ft suggests that the so-called bi-polar 
acoustic pulses generated at free metal-air interfaces 0] 
are particularly inefficient in magneto-acoustics since 
positive and negative parts in a bi-polar pulse cancel each 
other giving zero acoustic pulse area. Recently observed 
large-amplitude unipolar acoustic pulses with 3 ps dura- 
tion and amplitudes up to 1% generated in a thin cobalt 
transducer sandwiched between dielectric substrate and 
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FIG. 3: After application of large-amplitude 3 ps long acoustic 
pulse the magnetization vector initially prepared in state 2 
will decay (switch) into one of the neighboring energy minima, 
(a) a sequence of compressional pulses with r\ ac — —1.6%, 
induces the clock-wise rotation of magnetization vector: 2 — > 
4, 4 — > 2, etc. (b) a sequence of tensile pulses with r\ ac = 0.9% 
results into repetitive switching between two adjacent energy 
minima: 2 — > 1, 1 — > 2, etc. Dashed contours show the 
magnetization trajectories after excitation by a single pulse 
in the absence of damping. 



a layer of noble metal 2(J 21 1 are better suited for ex- 
perimental investigations in coherent magneto-acoustics. 

When using rectangular unipolar acoustic pulses the 
kick angle A9 ac is proportional to the product -q ac x r ac 
and thus can be increased by using larger strain ampli- 
tude r\ac or somewhat longer pulse duration r ac <C T prec . 

Figure 3 shows that both compressional (r) ac < 0, 
Fig. 3a) or tensile (rj ac > 0, Fig. 3b) rectangular unipolar 
acoustic pulses with r ac = 3 ps are capable of switch- 
ing the initial magnetization direction into the new min- 
imum, which represents the main result of this paper. 
For example, a three picosecond long compressional pulse 
with strain amplitude 1.6% switches the magnetization 
from minimum 2 into minimum 4 (2 — ¥ 4, see Fig. 3a). 
Application of a consequent identical compressional pulse 
bring the magnetization back into minimum 2 (4 — > 2, 
suggesting that a pulse train will result into clock-wise 
rotation of magnetization vector. In contrast, a sequence 
of 3 ps long tensile strain pulses will periodically switch 
the magnetization between minima 2 and 1 (2 — > 1 and 
1 ^ 2 see Fig. 3b). 

Therefore the results of our simulations suggest that 
a clean experimental demonstration of magneto-acoustic 
switching would necessarily imply a single-shot experi- 
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FIG. 4: Diagram of magneto-acoustic switching. Depending 
on the duration r ac and amplitude r\ ac of a unipolar tensile 
acoustic pulse, the magnetization vector initially prepared in 
state 2 undergoes the transition 2—^1 (dark blue area), 2—7-4 
(red area), 2 — ^ 3 (white area) or 2 — > 2 (no switching, light 
blue area). 



merit where the magnetic system is prepared in the same 
state before the action of the consequent acoustic pulse. 

The typical kick angle required for switching is about 
20° and the different threshold switching amplitudes and 
pathways for tensile and compressive pulses are caused 
by different hight of the potential barrier between the 
neighboring energy minima. The more general phase di- 
agram for magneto-acoustic switching is shown in Fig. 4, 
where the boundaries between different switching zones 
generally follow the l/r ac dependence, in agreement with 
the assumption that primarily the amplitude of out-of- 
plane acoustic kick A9 ac ~ rj ac x r ac = const determines 
the switching pathway. Similar analysis for the acoustic 
shear pulses leads to the same conclusions, in particular 
with respect to the acoustic pulse area and dependence of 
the switching amplitude on the acoustic pulse duration. 

It is quite remarkable that the threshold fluence for 
all-optical magnetization switching in a rare-earth com- 
pound GdFeCo induced by a single circularly polarized 
pump pulse almost does not depend on the optical pulse 
duration [12j. Longer optical excitation leads to longer 
acoustic pulses 22j, which are generated by thermo- 
elastic expansion of laser-heated GdFeCo. The analysis 
of thermo-elastic generation mechanism by longer optical 
pulses suggests that the absorbed laser fluence is propor- 
tional to the product of acoustic pulse amplitude by its 
duration, i.e to the acoustic pulse area. Therefore, the 
observed constant threshold fluence for all-optical mag- 
netization switching in GdFeCo could be explained by 
magneto-elastic mechanism. From the magneto-acoustic 
point of view a significantly lower magnetostriction co- 
efficient A 



(M)M S ~ 8 mT close to the compensation point [10j, thus 
giving the same order-of- magnitude prefactor A/M s in 
Eq. (|9]). The observed increase of switching fluence in 
GdFeCo with temperature [IH is consistent with the 
decrease of A/M s in rare-earth compounds for higher 
temperatures [l5|. Therefore, despite of being far away 
from quantitative modeling, all these arguments corrobo- 
rate the hypothesis [III that ultrashort acoustic phonon 
pulses may play an important role in the dynamics of 
all-optical magnetization switching in GdFeCo [13, EH]. 
Moreover, possible magneto-elastic contributions of heli- 
cally polarized acoustic shear pulses excited by circularly 
polarized electromagnetic pulses [3, [HI may solve the 
puzzle of an undefined long-living reservoir for angular 
momentum as highlig hted in the most recent systematic 
investigations 12j, |l3[ . 

To summarize, in this letter we have theoretically 
predicted the new mechanism of ultrafast non-thermal 
magneto-acoustic switching in Terfcnol-D. This phe- 
nomenon may open the door to ultrafast magnetic 
recording not relying on heating the magnetic material 
close the Curie point. 

Stimulating discussions with Stephane Andrieu and 
Karine Dumesnil and the financial support by Nouvelle 
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Loire are gratefully acknowledged. 



1CT 5 in GdFeCo [23| is balanced by a 



much lower room temperature saturation magnetization 
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